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Abstract 

The free energy of effective spin or "Polyakov line" models with a chemical potential, based on the U(iV) 
group, does not depend on the chemical potential. In a mean field-inspired expansion, we show how the 
condition of unit determinant, taking U(A0 to SU(/V), reintroduces the chemical potential, and allows us 
to express the free energy, as a function of mean field variational parameters, in terms of an expansion 
in the baryon (rather than the quark) fugacity at each lattice site. We solve the SU(3) mean field equations 
numerically to determine the phase diagram and compute observables. We also calculate the first corrections 
to the leading order mean field results, and find that these can significantly shift the endpoint of a line of 
first order transitions. The problem of deriving an effective spin model from full QCD is discussed. 
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I. INTRODUCTION 



Polyakov line or "effective spin" models, with lattice actions of the form 

d 

x k=l x 

are of interest as crude models of gauge theories in D = d + 1 dimensions at finite temperature 
and chemical potential [1]. Indeed, actions of this form can be extracted from QCD directly 
by integrating out most of the variables via a combined strong-coupling and hopping parameter 
expansion, while keeping the Polyakov line holonomies U x fixed, and therefore (1) is justified as 
an effective theory at least within the range of validity of these expansions. 1 At finite chemical 
potential /i the Polyakov line models have a sign problem, so that the usual Monte Carlo simulation 
is not directly applicable. There are, nonetheless, several different methods which can be used to 
solve this model. One of the earliest studies applied the complex Langevin equations to the SU(3) 
model [1-3]. A second method is the mean field approach, applied to the jU ^ case by Bilic et 
al. [2]. A third procedure, introduced in ref. [4], is to convert the partition function to a "flux" 
representation, which, in the SU(3) case, has been simulated via the Prokof'ev-Svistunov worm 
algorithm by Mercado and Gattringer [5]. Finally, the model can also be solved, at least in some 
parameter range, by the reweighting technique [6]. 

In this article we will revisit the mean field strategy, because there are certain aspects of that 
approach which we find illuminating. It is generally believed that the free energy of effective spin 
models based on the U(N) group do not depend on the chemical potential, for reasons discussed, 
e.g., in ref. [7]. We first rederive this result, in section II, in the framework of a mean field- 
inspired expansion. We then go on to show, in section III, how the restriction to a unit determinant, 
which converts U(A0 to SU(/V), not only reintroduces the chemical potential, but also converts the 
mean field formulation into an expansion in baryon fugacity. Numerical solutions of the mean 
field equations for the SU(3) are presented in section IV, and the phase diagram (projected to 
the /3 — /l plane) is obtained. We also display the effects of including the first correction to the 
mean field approximation. In section V we present some comments on the problem of extracting 
the appropriate effective spin model from full QCD, in the range of gauge couplings and quark 
masses of interest. Our conclusions are in section VI. 



II. U(A0 POLYAKOV LINE MODELS 

We will begin with models in which the effective spin (or "Polyakov line") variable U (x) is 
an element of the U(/V) group. As already noted, the chemical potential disappears from the free 
energy in this case, but the example will set the stage for the more interesting SU(N) models. 

Starting from the action (1), we mimic the mean field approach by first adding and subtracting 



Below we will refer to fi in eq. (1) as the "quark" chemical potential, while keeping in mind the fact that, in the 
hopping parameter expansion, jj. is actually related to the quark chemical potential of full QCD by a factor of 
inverse temperature. 
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constants u, v, which will eventually become variational parameters: 

S = p £ [(Trt/j - v + v) (TrU x+i -u + u) + (TrU x -u + u) (Trt/J^ - v + 1 

x,£ 

.V 

= -2/3JVMv + 2/3Jv£Trt/ A - + 2/3Jw£Trt/ A t + KrJ^Trt/, + e^Trt/J] +J . 

Here V is the lattice volume, <i is its dimensionality, and 

J = /3 £{ (TrC/t - v) (TrC/^ - «) + (TrU x - u) (Trt/J + - - v) } . 



We then have 



where 



S = -2j5dVuv + J^[A x TrU x + B x TtU^] + J 



A X = A = 2fidv + fce^ and B X = B = 2fidu + Ke~^ . 



(2) 



(3) 



(4) 



(5) 



Although A Xl B x are jc-independent constants, it is useful below to regard them as variables. This 
allows us to differentiate with respect to each of them, with the understanding that all the A X ,B X 
are set to A and B, respectively, after the differentiation. 

Ordinary mean field theory amounts to dropping J in the action and, in the absence of a chem- 
ical potential, setting, u = v = m, where m is the mean field. One then varies m to minimize the 
free energy. In our case, define 



Z mf = e- F "'f = e- 2 P dVuv H f dU x esp\A x TaJ x + B x TdJl] 



(6) 



and 



-AF 



Also defining the operator 



jDUe J expK^A./M/, + fl c Trf/*)] 
jDUcxp[£ x (A x TrU x + B x TrUl)} 



(7) 



J 

we have 



d d 

U, V, -r— , — — 

' ' dA' dB 



PL 

x,k 



eJlWTA'Tsl jDUexp[Z x (A x TrU x + B x TrU^)} 



dB, 



a a 



dA 



u + 



x+k 



dA, 



dB 



x+k 



(8) 



exp[-AF] 



fDUexp[Z x (A x TrU x + B x TrUl 
Next we need to evaluate the U(N) integral 

dU exp[AlW + STrC/ f ] , 



(9) 



\A X =A,B X =B 



(10) 
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which, by standard methods (cf. [8]), becomes an angular integration 

N N 



=1 2n N \ H - lH ]1 - JN 



exp 



A £ j+*+B £ g" 



m=l 



m=\ 



1 



N 

n 

A? 



2% 



exp 



/ 5 V 7 " / d V" /■ d(j) n 



This gives us 



IJUv [dBj J 2k 



^■j e h—tN e ji—jN 



exp 



2VA5 



z f=e^ 2pdVuv r\—e- e n 

■^m/ e 11 A7-i c «l— JJvll 



AH 



n=l 



<9A, J V dfl 



/o 



We now introduce rescaled variables 



u = € *V and v = e'V 
A x = (2/3Jv / + k)^=A^ 
5, = (2/3 du' + K)e-V = ti^T* 
d „ 3 



dA, 



3Ai 



dB x 



dB' x 



Then Z m y becomes 



/V 



Z f = e - 2 P dVLI ' v "n —£■ ■ £ 



-2PdVu'v' 



■m 
i 



n 

n=l 



(in-jn)p 



Ft e h-iN e h-iN ex P 



N 



Jn 



dB[ 



im=l m=l 



A / d \ 7 " / 5 ' 

x ii 



n=\ 



dB< 



k 



tyWx 



2^A 



At this point we note that, because of the £i l ...i N £j l ...j N term, 



N 



N 



Therefore 



1 N / 2 

Z f = e- 2pdVu ' v 'r\—£- £ TTl — 



n=l 



dA' 



dM 



J„ 



dB' x 



2jA' x B' x 



(16) 



As a function of the rescaled variational parameters u! ',v' ', Z m f is clearly, /i -independent, and of 
course it will remain /! -independent when F m f is minimized with respect to u', v'. Likewise, all /I 

dependence cancels in the J operator 



d d 



dB Y 



dA 



x+k 
d 



dA' 



u + 



u' + 



dA, 



— u 



dB 



x+k 



dA' x 



x+k 
d 



dB' , 

x+k 



(17) 



From this we can conclude that both F m f and AF, and therefore the free energy F = F m f + AF 
itself, are independent of the chemical potential /! in Polyakov line models based on the group 
U(A0. 2 

Before proceeding to SU(iV), we note that the expression for Z m f can be simplified a little 
further, using the identity 



d_d_ 
dAdB 

which is evident from the fact that 

/o 



/o 



2VAB 



/o 



2VAB 



2VAB 



Then, defining the derivative operator 



Dij(x) 



d§ ^e^+Be-^ 
2% 



(18) 



(19) 



(20) 



we may write 



Z mf = e- 2 ^ dVlN Y[^ fAy(x)/ [2v^] 



(21) 



and 



-AF 



1 



J\u'\ 



i _jL _jLi 

' dA' ' 3B> J Z, 



-mf 



mf 



(22) 



*a' x =a' \b' x =b> 

2 A slight subtlety is that at K = 0, the free energy depends not on uf, v' separately, but only on the product u'V = uv. 
Then one must appeal to the hermiticity of the action to set u = v. For any non-zero K and fi, however, there is no 
such degeneracy. 



Again, the jl -independence of the free energy is manifest. 



HI. SU(A0 POLYAKOV LINE MODELS 

We can convert the U(A0 models considered above to SU(A0 models by simply converting the 
U(A0 group integration in eq. (11) to an SU(A0 integration. To accomplish this, we only have 
to insert a periodic delta function into the angular integrations, which imposes the constraint that 

;1 = mod 2n. We use the identity 



N 

I 



1 



E ^ J = 2^ E ex P 



N 



E 0n 



(23) 



This introduces into each § n integration an additional factor of exp[is(j) n ]. Tracing through the steps 
of the previous section, we arrive at 



1 00 1 

7 -2pdVuvT~fJ_ y J_ r-r 

n= 1 \ 



2% 

x s=—°° 



( 3 



J \3B X J 1 (J<0 ) f^-V L V 



(24) 



Now expressing everything in terms of the rescaled variables of eq. (13), this becomes 



7 _ -IpdVu'v 1 FT J_ . _ . . J_ I y -sNu 1-T 
^mf — e LL N i tl i- l N t Ji-JN j7r < \Z^ e 11 



2% 



N 



s<0 n=l 



dA' x 



In 



S>0 n=\ 

in+U\ 



dA' 



S+jn 



BB' X 



dB' x 



(25) 



Defining 



Dij+ S (x) s>0 



(26) 



we can express Z m f compactly in the form 

CO 



X S=—c° 



(27) 



where we have also changed variables s — > — s in the sum. As before 



—AF 



j_/[«',v' ' 



- , mf 



mf 



(28) 



\A' X =A' fi' x =Bf 



This gives a formal expression for the full free energy, F(ji) = F m f(ji) + AF(jti) in terms of the 
variational parameters u', v', which should be chosen to minimize F(fi) . 

The mean field expression for the free energy F m f, as a function of the variational parameters 
u 1 ', v' (or equivalently A', B') has some features which are worth noting. In the first place, the mean 



This jU -independence was also demonstrated in the N = °° limit in [9]. 
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field partition function Z m t has now been expressed in terms of a product, at each site, of a fugacity 
expansion of the form 



£ e^det Dr s I [2VA>B i ] 



(29) 



Here we see that the quark chemical potential ji only occurs in the combination Njl, which is, in 
effect, the baryon chemical potential. So in fact we have an expansion in the baryon, rather than 
quark, fugacity. In ref. [10], the determinant in an expansion of this sort is referred as the "canon- 
ical determinant." The second point is that parameter s, originally introduced in the representation 
(23) of the periodic delta function, has now emerged as the baryon number (which, if negative, is 
the number of antibaryons) per site. 

Of course, one still has to minimize the free energy with respect to the variational parameters, 
and this will introduce some A/i-dependence into the canonical determinants. Strictly speaking, it 
is the mean field expression of the partition function as a function of (freely varying) parameters 
u', v' which has the form of a fugacity expansion. 

Successivejmprovements to the leading mean field result would be obtained by expanding the 
operator exp[/] in a Taylor series. In the case that K = 0, and /3 is so small that the minimum 
free energy is obtained at u' = v' = 0, then the Taylor series simply generates the strong-coupling 
expansion. At larger /3 and K, the series also generates corrections to the leading mean- field result. 
We will compute the effect of the leading correction in the next section. 

At this point, we should draw attention to the similarities and differences between our approach 
and the much earlier work of Bilic et al. [2]. The starting point of the mean field treatment in [2] 
was the action (4) without the /-term. The SU(3) group integral was expanded as a power series 
in A, B, and for this reason it was not obvious that the partition function is an expansion in baryon 
fugacity, arising from the unit determinant condition. In the next section we determine the phase 
diagram (for both real and imaginary /i), which was not displayed in [2], and work out leading 
corrections to the mean field result. 



IV. NUMERICAL RESULTS FOR THE SU(3) POLYAKOV LINE MODEL 



We will now specialize to SU(3). From eq. (27), we see that the mean field free energy per 
lattice site at ./V = 3 is 



fmf = 2/3dwV-log 



£ e 3 ^det D^/opv^ 7 ] 



(30) 



where we have dropped an irrelevant constant. In numerical work we cannot sum s over the full 
range [— °°,°°], so it is necessary to cut off the sum at some maximum baryon/antibaryon number 
Smax per site 



f mf ^2pdu'v'-\og[G(A',B')} 



s max 



G{A\B')= £ e^det Djf I^iVa^B 1 



(3D 



and of course it is important, when computing observables, to check sensitivity to the cutoff. We 
will return to this issue below. 
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Minimizing the free energy with respect to the variational parameters u' ,V, or, equivalently, 
with respect to A' = 2/3 dv' + K,B' = 2f3du' + K, leads to two equations 



B'-K 



1 



dG 



2/3 d G(A',B')dA' 
A'-K 1 dG 



2/3 d G(A',B')dB' 








(32) 



whose roots may be determined numerically. 4 At the minimum, we can regard A' = A'(/3, K, [l) 
and B' = 5'(/3, as functions of the parameters of the theory. 

Apart from the free energy itself, the observables of interest are Tr[t/],Tr[L^], and the baryon 
number density n (baryon number per lattice site). The latter is given by 



(n) 



dfmf 

dm 

s max i- 

— — — T se 3s ^det D^IollVAW] 
G(A>,B>) s _ u s L 3 J 

1 fdA'(p.,K,n) d dB'(B.,K,n) d 



M + dH dff) fmf{A ' B) 



(33) 



where it is understood that the derivative is taken at the point where f m f(A r ,B r ) is minimized. But 
at this point, the first derivatives of f m j with respect to A' and B' vanish. Therefore 



{n) = g(A5o se3s " det r^ 7 « t 2 ^ 7 ] 

V * / S — S run v 



(34) 



From (6) we see that 



= A log G{A', B') 

= e-^\o g G(A>,B>). 

At the minimum of the free energy, determined by the roots of (32), this simply becomes 

1 



(Trf/> = -£(Trt/,) = e -'V = u 



(35) 



(36) 



and likewise 



( TrC/t ) = ^L( Tr ^) = ^ 



V = V . 



(37) 



This is, of course, reminiscent of the standard mean field approach to a spin system, in which the 
variational parameter becomes the average spin. It must be understood, however, that due to the 



4 We have found Mathematica convenient for this purpose. 
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complex weight there is no constraint that the "average" values of TrU and TrU^ are necessarily 
bounded by Trl. 



k=0.02, (i=1.2 



k=0.02, (1=1.2 



B <Tr[U + ]> 
.6 

4 - 

2 - 

1 - 
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(a) Polyakov lines vs. j3. 



0.1 0.11 0.12 0.13 0.14 0.15 0.16 
P 

(b) - free energy/site vs. j3. 



FIG. 1. Observables vs. j8 at fixed K = 0.02 and /I = 1.2, evaluated at mean field level for SU(3). 



We now have all the tools needed to evaluate observables and map out the phase diagram. 
Figure 1 shows a typical result for (Trt/), (Trt/^) and the mean field free energy per site f m f, as a 
function of /3, at fixed K = 0.02 and chemical potential /! = 1.2. There is a clear first order phase 
transition at /3 = 0. 1257. As the chemical potential is increased at fixed K = 0.2, the discontinuity 
at the transition decreases, until it disappears altogether at /! = 1.67. At larger /i, there is only a 
crossover. 

Repeating this procedure, we can map out the region of first order transitions in the /3,/i, K 
parameter space. In Fig. 2 we show sample first-order phase transition lines in the /3 — /i plane at 
K = 0,0.02,0.04,0.05,0.059. At K = the transition, at /3 = 0.1339, is of course independent of 
fi. At fixed, finite K the transition line terminates at some value of jl, and this termination point 
happens at smaller and smaller values of jLi as K increases. The transition line shrinks to a point at 
/l = for K = 0.059, and beyond this value of K there are no further transitions. 

We can also solve the mean field equations for imaginary /i. The results for several values of fc 
are shown in Fig. 3. The continuity of first order transition lines, as /i varies from real to imaginary 
values, ties in with the considerations of ref. [11]. 

Figure 2 can be compared directly to the phase diagram recently obtained by Mercado and 
Gattringer [5] via a Monte Carlo simulation in the flux representation. The two diagrams are 
qualitatively, and even quantitatively, very similar. The main difference is that we only show first 
order transition points, and most of these are found in ref. [5] to be crossover points, rather than 
first order transitions. According to Mercado and Gattringer [5], the endpoint of a line of first 
order transitions, at a given K, occurs at a much smaller value of /! than we find in our mean field 
calculations. So an interesting question is whether inclusion of higher order corrections, beyond 
the leading order mean field result, would bring our endpoint s to smaller values of /i, in closer 
agreement with [5]. We will turn to this question in subsection IV B below. 



A. Effect of the baryon number cutoff 

The data displayed above was obtained using a cutoff s max = 4 in the sum over baryon number, 
but the results shown are quite insensitive to increasing the cutoff to s max = 6, and even to decreas- 
ing the limit to s max = 2. The reason for this insensitivity is that the phase transitions occur at 
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0.14 
0.135 
0.13 
0.125 
CQ_ 0.12 
0.115 
0.11 
0.105 
0.1 

0.5 1 1.5 2 2.5 

FIG. 2. Phase diagram of the Polyakov line model (1) for the SU(3) group, obtained via mean field methods, 
in the j8 — fi plane at several values of K. The lines indicate first order transitions. Beyond K = 0.059, there 
are no transitions at any value of }x. 

values of the baryon number density which are very small compared to the cutoff. Only when the 
chemical potential is raised to values such that the number density becomes comparable to s max 
does the cutoff dependence become apparent. To illustrate this dependence, we fix /3 = 0.1257 
and K = 0.02 (where we have found a transition at /i = 1.2), and compute the Polyakov lines and 
number density over a wider range of /i. 

The results, for [i < 10 and s max = 2,4, are shown in Fig. 4. We see that (TrU) and (TrU^) are 
comparable to one another and of O(l) until (n) approaches the cutoff in s. Beyond that point, 
(TrU) falls exponentially as and (TrU^) diverges as e^, exactly as in the U(/V) theory, and 
the results are no longer valid for the SU(A0 case. When (n) saturates the cutoff then, in order to 
probe a larger range of /i, it is necessary to increase s max . For the purpose of determining the phase 
diagram, however, s max = 4, which can be interpreted as a limit of no more than four baryons per 
lattice site, appears to be more than sufficient. 

B. The leading correction to the mean field free energy 

Going back to eq. (7), we have 




The product is over all links, where 

J x , k = /3 { (Trt£ - v) (TrU x+i -u) + (TrU x - u) (Trt/J + - - v) } , (39) 



1 1 1 1 


F -— ^ 
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Transition points for imaginary u 
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FIG. 3. Some transition points for Polyakov line models in the /$ — Im^) plane, for imaginary values of the 
chemical potential, at several values of K. Filled circles indicate first order transition points, open circles 
indicate a crossover. 



and the () m y notation denotes the expectation value with respect to the mean field action, as in (7). 
The expansion of exp[7] generates products of terms such as Ji x Ji 2 ...Ji n , where the /, denote links, 
some of which may be the same. Because (Tr[U]) m f = u and (Tr[t/']) m y = v, it is clear that the 
expectation values of such products are only non-zero if each endpoint of a link /, appearing in 
the product is also an endpoint of at least one other link appearing in the product. The simplest 
product whose expectation value is non-vanishing, containing the minimum number of 7 factors, 
is simply the product of 7/7/ on the same link. Therefore, to leading order, we approximate 




«n(i+^*w). (40) 



ii 



k=0.02, B=0.1257 
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10 



0.1 



0.01 



s max~^ 

Smay = 4 



1000 



100 



+ 
3 



2 4 6 

(a) (Trt/) vs. jU. 
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(b) (Trt/ 1 ") vs. jU. 
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(c) number density vs. jx. 

FIG. 4. Observables (Tr{7), (Tr£/ f ) and (n) vs. /I at fixed K = 0.02 and j3 = 0.1257, for two vafues of the 
cutoff s max . Note that these observables are independent of the baryon/site cutoff s max , until a little beyond 
jU = 4, which is well past the value of jU at the first order transition. 



Now 



+ (1W; - v) (Trt/, - u) (TrU x+l - u) (TrC/J + . - v) + /i.e. 



= 2/3 : 



(Tr£/Trt/) m/ -a 2 ) ((Trt^Trt/^-v 2 ) + ((Trt/Trt/ t ) m/ - uv 



and we use 



2/3 2 [(((^Trt/) 2 ) m/ - M /2 ) ((( e -M W t)2) 



+ f ((^Trf/)( e -^Trt/ t )),„ / - wV 



1 



/2 



(^Trf/) m (e-^Trt/ t )' 



«/ G(A',fl') V^A'y \dB' 
Putting all the pieces together, the free energy per unit volume is 

f(A',B') = 2pduv-f(A',B'), 



G(A',B') . 



(41) 



(42) 



(43) 
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where 



f = logG(A',B')+dlog 



1 d 2 G 
GdA' 2 



,12 



1 d 2 G 
GdB' 2 



12 



+ 



1 d 2 G 
G dA'dB' 



u'v' 



(44) 



and G(A',B') is as defined in (31). Note that the terms inside the logarithm, which correct the 
leading mean field expression, depend on fluctuations around the mean field values. 

The variational parameters A',B' are again derived by minimizing f(A',B'), which implies 



B' 



K 



2fid 
A'-K 



dA 
d 



2pd dB 



7/ = ° 
7?=° 



(45) 



whose roots may again be determined numerically. It is also still true that u = (Trf7),v = (Trt/^), 
which can be seen as follows: Define 



f DUe J exp £(ATr£/ x + flTrf/J 



Then it is clear that 



(46) 



(47) 



Applying the first of eqs. (45), and the definitions (13), we arrive at u = (Trf7). In the same way, 
we can show that v = (TtW). Thus the correspondence between the variational parameters u,v 
and the observables (Trt/), (TtU^) is maintained exactly, in fact to all orders beyond the leading 
mean field expressions. 

We can now study how inclusion of the leading correction will modify the phase diagram 
shown in Fig. 2. It turns out that the location of the phase transition points changes very little. 
Generally, at fixed k, /i, the value of /3 at the transition changes by less than one percent. What 
does change significantly are the endpoints of the first-order transition lines. For example, at 
K = 0.02, the endpoint of the transition line was at /i = 1 .67, /3 = 0.1213. Inclusion of the first 
correction brings the endpoint down to /i = 1.38, /3 = 0.1249. The free energy at lowest order 
(mf), and the free energy after inclusion of the first correction (mf+nlo) is shown in figure 5. 
The free energy changes substantially, but the transition point hardly at all (from /3 = 0.1243 to 
/3 = 0.1249). However, at /l = 1.38, the order of the transition changes, from first order in the 
leading mean-field approximation, to a sharp crossover when the first correction is included. In 
figure 6 we show a closeup of the (Trf/) in the near neighborhood of the transition in both cases. 

We also find that at k = 0.04, the endpoint of the line of first order transitions moves from 
ji = 0.87, /3 = 0.121 1 to n = 0.46, /3 = 0.1246. At K = 0.045 the line of transitions shrinks to a 
point, at ji = 0, /3 = 0.1245. Beyond K = 0.45, there are no transitions. 

So the first correction to mean field is taking us in the right direction, in the sense of bringing 
the endpoint of the first order transition line to smaller values of jU. Mercado and Gattringer [5] 
find that the endpoints of the first-order transition lines are located at yet smaller values of /i. It 
would be interesting to see if the next higher-order corrections generated by exp [J] would bring the 
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FIG. 5. — F in lowest order mean field theory, and in mean field + next leading order. Inclusion of the next 
leading order can change a first-order transition to a crossover, as seen clearly in Fig. 6. 
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FIG. 6. Closeup of (TrU) in the transition region, at fc = 0.02 and pL = 1.38, showing the effect of inclusion 
of the leading correction to the mean field free energy. With the inclusion of the leading correction, this 
is the endpoint of the k = 0.02 line of transitions, down from the value jU = 1.67, which is the endpoint 
without the first correction. 



endpoints still closer to the endpoints found in ref. [5]. We leave this exercise for a future study. 



V. EFFECTIVE SPIN MODELS AND FULL QCD 

It seems to be easier to solve effective spin models at finite chemical potential, by a variety 
of methods, than to solve full QCD at finite chemical potential. This means that if we knew the 
effective spin models corresponding to full QCD at relevant points in the plane of temperature and 
quark chemical potential, then by solving the effective models we could determine the QCD phase 
diagram. We know how to derive the effective spin model in the strong coupling and hopping 
parameter expansions; for jU = this has been done in [12, 13], and for /! ^ in [6]. The problem 
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is to go beyond these expansions, to weaker couplings and light quark masses, at /i 7^ 0. 5 

In principle, the effective Polyakov line model is derived from full QCD by integrating out 
the quark and gauge field variables, under the constraint that the Polyakov lines are fixed. It is 
convenient to impose a temporal gauge on the periodic lattice, in which all timelike links are set 
to the unit matrix except on a single time slice, t = say. Then the effective theory, at chemical 
potential /1 = 0, is defined by integrating over all quark fields and link variables with the exception 
of the timelike links at t = 0, i.e. 

Z(p,T,m f ) = JdU (x,0) J DU k DyDy e s Q CD 

= J DU (x, 0) e Se ff [u °' u ° ] , (48) 

where /3 is the gauge coupling, T = l/N t is the temperature in lattice units with N t the lattice 
extension in the time direction, and nif represents the set of quark masses. Because temporal 
gauge has a residual symmetry under time-independent gauge transformations, it follow that S e ff 
is invariant under Uq(x,0) — » g(x)i7o(x,0)g^(x), and therefore can depend on the timelike links 
only through their eigenvalues. This just means that S e ff is a Polyakov line action of some kind. 

Let Sq CD denote the QCD action with a chemical potential, which can be obtained from Sqcd 
by the following replacement of timelike links at t = 0: 

S I1 QC d = s QCd[Uo(x,0) ^e^t/ (x,0),t/ t (x,0) e'^U^O)] . (49) 

The effective Polykov line action, at finite chemical potential is defined via 

Z(n,p,T,m f ) = JdU (x,0) J DUkDyfDyr e s Qcv 

= J DU Q (x,0) e S eff [Uo > u ° ] . (50) 

The integration over [4, y7, i/a can be carried out in a strong gauge-coupling and hopping pa- 
rameter expansion, to obtain S e ff and It is not hard to see that each contribution to S e ff in the 
strong coupling + hopping parameter expansion of (48) maps into a corresponding contribution to 
Seff, in the expansion of (50), by the replacement 

Ux ^ e NtH Ux 5 uZ^e-^ul, (51) 

where we have identified U x = £/q(x,0). Since this mapping holds to all orders in the strong- 
coupling + hopping expansion, it is reasonable to assume that it holds in general, i.e. 

Suffix, Ul] = S eff [U x e N '»U x , Ul e-^Ul) . (52) 

Equation (52) is a rather trivial, but potentially powerful identity. It suggests that if, by some 
means, one could obtain S e ff at fixed {[5,mf, T} and chemical potential ji = 0, then one would 
immediately also have the effective action S^ f at the same set of parameters {/3,m/, T}, but any 
chemical potential /1, by the replacement shown. Unfortunately, there is some degree of ambiguity 
in S e ff at /! = 0. Suppose we have some ansatz for S e ff[U x , UJ], depending on some small set of 
parameters, which we would like to fix by comparing to the full theory at ji = 0. The problem 



For efforts at deriving the effective Polyakov line model in pure gauge theories, c.f. [14] and references therein. 
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is that whatever ansatz we make for S e ff, there is another form which is identical to that ansatz 
at /i = 0, but differs under (51). In the case of SU(2), for example, we could use the fact that 
TrU^ = TrU to replace Tr£/J by Trt/ X (or vice versa) in S e ff. Obviously, this will not affect 
anything at /i = 0, but will result in a different theory at non-zero /i. Similarly, in SU(3), the 
identity 



+ 1 

x 2 



(TrC/ x ) 2 -Trt/, 2 



(53) 



allows us to replace Trt/J" everywhere in S e ff by the right hand side of (53), but this again produces 
quite a different theory at ji ^ under the rule (52). Of course a similar identity holds for Trt/ X , so 
we can convert the original S e ff to another theory which may be symmetric in U x , U], but which 
has quite a different extension to finite chemical potential. 

It may be possible to overcome this ambiguity, however. Suppose we take the timelike link 
variables at t = to be U(3), rather than SU(3) matrices. Then the ambiguity due to (53) is no 
longer present, but the effective spin theory still only depends on the eigenvalues of the U(3) ma- 
trices. Then let us suppose that we have some reasonable ansatz for S e ff in a physically interesting 
range of parameters p,m,f,T, e.g. 

S eff = £7(x-y)Tr[£/jTr[t/;] + £/(x-y) (Tt[U x ]Tt[U y ) +Tr[^]Tr[t//]) + £V(tf„E£) , 

x,y x,y x 

(54) 

where /(*), J'(x) are parametrized by a few constants (such as nearest and next -nearest neighbor 
couplings), and V(U X , U%) can be limited to a few terms involving the characters of U(3). In that 
case, the effective spin model is specified by a handful of constants {cj}, which of course depend 
on {j8,m/,r}. 

Since there is no sign problem at /! = and U x = Uo(x,t = 0) e U (3), it should be possible to 
numerically simulate both the effective theory and the full theory. Then one can imagine a number 
of strategies for obtaining the constants {cj}. One possibility is to simply calculate an appropriate 
set of observables in both theories (Polyakov lines in various representations and Polyakov line 
correlators), and fix the set of constants {cj} in S e ff so that the two theories yield the same re- 
sults. Or perhaps some variant of the inverse Monte Carlo method could be applied [14]. A third 
procedure is inspired by a recent study of the Yang-Mills vacuum wavefunctional [15]. The idea 
is to select a finite set of M timelike link configurations 

{U® = U®(x,t = 0) e U(3),i=l,2,...,M} , (55) 

where each member of the set specifies the timelike link variables at every spatial site x and 
t = 0. Then the Monte Carlo simulation of the full theory proceeds in the usual way, except that 
on the t = timeslice, one member of the given set of timelike link configurations is selected by 
the Metropolis algorithm, and all timelike links on that timeslice are updated simultaneously. Let 
be the number of times that the ?'-th configuration is selected by the algorithm, and N tot = M - 
Then it is not hard to show that 



exp 


'Seff[UU[ 


exp 


's eff [uu)[ 



lim . (56) 
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Information derived from a number of such simulations, each using a different set of configua- 
tions at t = 0, can in principle completely determine the {c/}. However, since the {Ni} vary 
exponentially with S e tf, the variation of S e ff within a given set must be kept relatively small, i.e. 
8S e ff ~ 5 — 7, in order to ensure a reasonable acceptance rate for all members of the set. For de- 
tails of the algorithm, and its application to the vacuum wavefunctional of pure Yang-Mills theory, 
cf. [15]. 

Once the set of constants {cj} is found, by whatever method, the effective theory at finite 
chemical potential, S^//' ^ or an y ^ but *- ne same set {P, m f, T}, is given by the identity (52). The 
final step is simply to note that SU(3) C U(3), so that the theory we want, is obtained by 
restricting the U x matrices to the SU(3) subset. Equivalently, since we can always express the U(3) 
matrices as 6 

U x = e xp[iO x }U° U (3) , U^=cxp[-ie x ](uf {3) y (57) 

the conversion from S e ff to is obtained by setting X = —iN t jl. 

With the effective Polyakov line model S^y in hand, the theory can be solved by the mean 
field approach discussed above, or by other methods such as complex Langevin [1-3], the flux 
representation + worm algorithm [5], or reweighting [6]. To check that the method is working at 
jli 7^0, one would compare full QCD with the effective spin model at, e.g., small or imaginary /i, 
where the /l -dependence of the full theory can be obtained by other means. 

This approach can be expected to break down at sufficiently large /I. At some point, terms in 
the potential involving high powers of U x and U^, which might be negligible for computing ob- 
servables at /i = because they are multiplied by very small coefficients, could become important 
under the replacement (52). To what extent this effect will inhibit the study of the phase diagram 
remains to be seen. 

There is no doubt that determining the set of constants {c,} in full QCD would be computa- 
tionally demanding. As a first step, it may be worth trying to extract the effective spin theory from 
gauge theories with scalar, rather than fermionic, matter fields. 

VI. CONCLUSIONS 

The mean field expansion for effective spin models with a chemical potential turns out to have 
an interesting structure. The constraint taking U(A0 to SU(/V) is responsible for the /i-dependence 
of the free energy, and this constraint introduces an infinite sum whose index, as it turns out, can 
be interpreted as the baryon number at each site. The partition function can then be formally 
expressed in terms of a baryon fugacity expansion. 

If we ignore the distinction between first-order and crossover points, then even the lowest order 
mean field equations do a reasonably good job of accounting for phase structure. The main error 
lies in the location of the endpoints of first-order transition lines, which occur, for fixed K, at 
higher values of /i than those determined by other methods. The first correction to the mean field 
result moves those endpoints in the right direction, i.e. to lower values of /I. It remains to be seen 
whether realistic results for the endpoints would be obtained from still higher orders in the mean 
field expansion. 

We have also commented on the problem of deriving effective spin models from full QCD, 
and on certain subtleties associated with continuing those models from zero to finite chemical 

6 Allowing for the Z3 subgroup of SU(3), the angle 9 X can be restricted to the range [0,2tt/3). 
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potential. We have suggested that a method which was previously applied to study the Yang-Mills 
vacuum wavefunctional may be useful in this context, and hope to discuss this further at a later 
time. 
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